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Abstract 

I give a brief outline of the theoretical framework for the modern treatment of the 
strong interaction effects in heavy quark decays, based on first principles of QCD. 
This mo del- independent approach is required to meet the precision of current and 
future experiments. Applications to a few problems of particular practical interest 
are reviewed, including the precise determination of V c b and V u b- I emphasize the 
peculiarities of simultaneously accounting for the perturbative and power-suppressed 
effects necessary for accurate predictions. 
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1 Introduction 



A key role in exploring the Standard Model is played by studying electroweak heavy 
flavor decays. It was realized 20 years ago that the strong interaction effects in 
heavy flavor hadrons can be treated within QCD. Yet the full power of theoretical 
methods acquired in QCD was applied here only recently. They were developed 
along two main directions, 'symmetry-based' and 'dynamical'. These two lines in 
the heavy quark theory were the counterparts of the basic theoretical strategy in 
studying strong interactions: isotopic invariance and chiral symmetry on the one 
hand, and asymptotic freedom on the other. 

In heavy quark physics, the early period to the end of the 80s saw mostly the 
dynamical approach applied at a simplified 'intuitive' level. The nonperturbative 
effects were often thought to be small even in the decays of charm particles. The 
following few years were dominated by 'symmetry' considerations; the operating 
language for those analyses was the so-called Heavy Quark Effective Theory (HQET), 
which incorporated some basic elements of the general heavy quark expansion in 
QCD (HQE) but was limited only to certain classes of processes. 

Finally, over the last few years a consistent well-defined dynamical approach 
has been developed, which automatically respects the heavy quark symmetries in a 
manifest way. Here the most precise determinations of \V c b\ and \V u j,\ were made. 

The main effects in a weak decay of heavy quarks Q originate from distances 
~ 1/rriQ <C 1/Aq CD . Since Q s (mg) <C 1 they are tractable through perturbation 
theory. The QCD interaction becomes strong only when the momentum transfer is 
much smaller than the heavy quark mass, k <C rriQ. Two basic ingredients of HQE 
are thus elucidated: 

• The nonrelativistic expansion, which yields the effects of 'soft' physics in the form 
of a power series in I/ijiq. 

• The treatment of the strong interaction domain based on the Operator Product 
Expansion (OPE). 

Unless an analytic solution of QCD is at hand, these two elements appear to be 
indispensable for heavy quark theory. 

The general idea of separating the two domains and applying different theoretical 
tools to them was formulated long ago by K. Wilson in the context of problems in 
statistical mechanics; in the modern language, applied to QCD it is similar to lattice 
gauge theories. The novel feature we face in the theoretical analysis of beauty de- 
cays is that they often allow - and even demand by virtue of the existence of precise 
experimental measurements - rather accurate predictions, requiring a simultaneous 
treatment of perturbative and nonperturbative QCD effects with enough precision in 
both. This problem is not new; the theoretical framework has been elaborated more 
than 10 years ago 0, but its phenomenological implementation was not mandatory 
until recently. Failure to incorporate it properly leads to certain theoretical para- 
doxes and, unfortunately, some superficial controversy in the numerical estimates in 
the literature. 
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Figure 1: Exclusive B — > (a) and generic (b) semileptonic decays. 

With significant progress made over the last years, the theory of the heavy flavors 
is still not a completed field and is undergoing to further extensive development. I 
will focus on a few selected topics that illustrate the theoretical framework, and 
review the overall status of the heavy quark expansion, with the main emphasis on 
the qualitative features. Some important theoretical applications are presented in 
the lectures by C. Sachrajda |J (exclusive decays), A. Ali || (rare b decays) and 
M. Gronau || (CP violation) (these Proceedings). Additional theoretical aspects 
are covered in the summarizing contribution by G. Martinelli ||, where recent ex- 
perimental data are also discussed. 

2 Semileptonic decays 

The Q CD-based heavy quark expansion can equally be applied to all types of heavy 
flavor transitions. Semileptonic decays are the simplest case and I shall devote most 
of the attention to them; for practical reasons I focus on b — > c transitions. A brief 
discussion of the b — > u decays will be given later. 

A typical semileptonic decay is schematically shown in Figs. 1. Generally, two 
types of decay rates can be singled out: inclusive widths where any combination of 
hadrons is allowed in the final state, and exclusive decays, when a transition into a 
particular charmed hadron is considered, usually D or D*. 

2.1 Inclusive semileptonic width 

The semileptonic width of a heavy quark has the form 



where zq is the known phase space suppression factor and ae generically includes all 
QCD corrections. In the heavy quark limit the difference resulting from using the 
quark mass m& and the meson mass Mb in eq. ([!]) disappears: 




(1) 



(M B - m h )/m b ~ l/m b 



(2) 



For the actual b quark M%/m\ amounts to a factor of 1.5-2, which formally con- 
stitutes a power-suppressed effect. This demonstrates the necessity of a systematic 
control of nonperturbative corrections even in decays of beauty particles. 

The central result obtained by direct application of OPE to the inclusive decay 
widths in QCD is the absence of 1/mg corrections [0] - in contrast with the presence 
of such terms in the hadron masses. The physical reason behind this fact is the 
conservation of the color flow in QCD, which leads to the cancellation of the effects 
of the color charge (Coulomb) interaction in the initial and final states. In terms of 
nonrelativistic quantum mechanics (QM), it is the cancellation between the phase 
space suppression caused by the Coulomb binding energy in the initial state, and the 
Coulomb distortion of the final state quark wavef unctions. The inclusive nature of 
the total widths ensures that they are sensitive only to the interaction on the time 
scale ~ 1/AE ~ 1/rrib. The final-state-interaction effect is thus not determined by 
the actual behavior of the strong forces at large distances, but only by the potential 
in the close vicinity of the heavy quark. The cancellation therefore occurs universally, 
whether or not a nonrelativistic QM description is applicable. 

The leading power corrections start with terms 1/m 2 ; they were calculated in [7] 
and are expressed in terms of the expectation values of two operators of dimension 
5, which have a transparent QM interpretation: 

/4 * ^(B\f-a^b\B) «-> (B\a b ■ g s H g \B) ~ ^(M 2 B , - M B ) ~ 0.35 GeV 2 (3) 

/4 * ^{BWDf^B) <-> (B\f\B) . (4) 

The value of /i 2 is not yet known directly; a model-independent lower bound was 
established in || ||: /x 2 > /jlq ; this puts an essential constraint on its possible 
values. This bound is in agreement with QCD sum rule calculations JTU], yielding 
a value of about 0.5 GeV 2 and with a more phenomenological estimate [Til]. In the 



absence of gluon corrections, as in simple QM models, the expectation value /i 2 
would coincide with the HQET parameter —Ax; they are different, however, in the 
actual field theory, where both /i 2 and \j? g depend on the normalization point. 

Including the nonperturbative corrections, the semileptonic width has the follow- 
ing form % [TJ, @, 0] 



GWb \V cb ? L(l- - 2 ( 1 - 4 V 4 - ?-4 1} + ...}■ (5) 



1927r3 rc01 | V 2m 2 ) \ ml ml 3 vr 



The 1/m 2 corrections to T s \ are rather small, about —5%, and increase the value of 
\V c b\ by 2.5%; the impact of the higher order power corrections is negligible. 

Good control of the QCD effects in the inclusive semileptonic widths provides 
the most accurate direct way to determine \V c b\ in a truly model-independent way. 
It sometimes faces a traditional scepticism: which numerical value must be used 
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Figure 2: Perturbative diagrams leading to the IR renormalon uncertainty in niQ° e 
of the order of Aqcd- The contribution of the gluon momenta below itlq expresses 
the classical Coulomb self-energy of the colored particle. The number of bubble 
insertions into the gluon propagator can be arbitrary. 



for nit, and m c ? This practical problem has deep roots; failure to understand them 
is the major source of controversy about masses and inclusive widths found in the 
literature. It will be briefly discussed below. In reality, the precise value of m& is 
not too important, since the b — > c width depends to a large extent on the difference 
nib — m c rather than on m& itself; the former is constrained in the HQE: 

M B + 2>Mb* M d + 3M d * 2 ( 1 1 \ 
m b -m c = i i + /H^-^J+-- 3 - 50GeV - ( 6 ) 

It also independently enters lepton spectra in semileptonic decays ]14| and can be 
extracted from the data |L5] . Numerically [§, [L6| , a change in m b by 50 MeV leads 
only to a 1% shift in \V c b\. 

Heavy quark masses 

The controversy about m& is due to the fact that HQET was popularly based on 
the so-called 'pole' mass of the heavy quarks. Not only was it a starting parameter 
of the HQET-based expansions, it is this pole mass that one always attempted to 
extract from the experimental data. It turns out, however, that the pole mass of the 
heavy quark is not a direct observable and its definition suffers from an irreducible 
intrinsic theoretical uncertainty of order Aqcd FJ - 



At first sight this looks paradoxical and counter-intuitive: for example, the value 
of m e quoted in the tables of physical constants is just the pole mass of the electron. 
In QCD there is no 'free heavy quark' particle in the physical spectrum, and its pole 
mass is not well defined. The problems facing the possibilities to extract the pole 
mass from typical measurements were illustrated in Refs. []T8|| and JOI . 

The physical origin of the uncertainty SniQ° le ~ Aqcd is the gluon Coulomb self- 
energy of the static colored particle. The energy stored in the chromoelectric field 
inside a sphere of radius R 1/ nig is given by 



a s (R) 1 

Hi r a x cx cons i — 

'l/mi,^\x\<R 

The pole mass assumes that all energy is counted, i.e. R — > oo. Since in QCD the 



SE Co niomb{R) oc / d?x oc const - ^_^_}_ _ (7) 

Jl/m h ~\x\<R 7T R 
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interaction becomes strong at R ~ 1/Aqcd, the domain outside R would yield an 
uncontrollable and physically senseless contribution to the mass ~ Aq CD |L7| . 



Being a classical effect originating at a momentum scale well below mg, this 
uncertainty can be traced in the usual perturbation theory, where it manifests itself 
in higher orders as a so-called 1/rriQ infrared (IR) renormalon singularity in the 
perturbative series for the pole mass p0| , see Fig. 2. 



Nonetheless, the inclusive widths can be theoretically calculated since they are 
governed, instead, by well-defined short- distance running masses mg(/j) with the 
Coulomb energy originating from distances > peeled off. It is precisely this 
short-distance running mass that can be extracted from experiment with, in prin- 
ciple, unlimited accuracy: the pole mass does not enter any genuine short- distance 



observable at the level of nonperturbative corrections |f20 . 

Applied to the inclusive widths, it suggests certain information about the impor- 
tance of higher order perturbative corrections: if masses entering eq. are the pole 
masses, the perturbative series 

rP ert = r ae pcrt = r„ (l + ^(ajn) + a 2 {a s /-nf + ...) (8) 

is poorly behaved, with coefficients factorially growing, which makes the radiative 
correction factor uncalculable in principle with an accuracy ~ Aqcd/^6- In contrast, 
if one uses the short-distance masses, the higher-order corrections become smaller 



and the factor ae pert becomes calculable with the necessary precision p0| ^2 . 

This seemingly academic observation, in reality proved to underlie the pattern of 
the corrections from the very first terms. Remarkably, the actual model-independent 
calculations of T s \ through observables measured in experiment are very stable 
against perturbative corrections. Including 0(a 2 s ) terms in the extraction of the 



b pole mass from, say, the e + e~ — > bb threshold region p3| noticeably increases 
its value. However, the parallel perturbative improvement in calculating the width 
yields an essential suppression of the perturbative factor, 33, so that the two effects 
offset each other almost completely ||16| . 



This conspiracy is not unexpected: the appearance of large corrections at both 
stages is an artefact of using the ill-defined pole mass in the intermediate calculations. 
The situation is peculiar since the actual nonperturbative effects appear only at the 
level l/rriq, whereas the pole mass is infrared ill-defined already at an accuracy of 



1/rriQ. The failure to realize this fact led to the superficial suggestion |23] that even 
in beauty particles the perturbative corrections may go out of theoretical control; a 
more careful analysis ^] showed that this is not the case. 

Moreover, the OPE requires using short- distance running masses m;, jC (/i) normal- 
ized at fi ~ 1 GeV |20| . It has been done in |16| and demonstrated that neither these 



masses nor the perturbative corrections to the width as show significant contributions 
from higher orders. 

To summarize, the idea that the perturbative corrections in the extraction of | V c b\ 
from T s i(B) are large comes from an inconsistent usage of ill-defined pole masses: 
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• It is 'difficult to extract' accurately ml° le from experiment; in any given calcu- 
lation it is easy to identify the effects that were left out, which can change its value 
by ~ 200 MeV. This uncertainty leads to a 'theoretical error' 8j in T a \{B) of ~ 10%. 

• When routinely calculating r s i(5) in terms of the pole masses, there are sig- 
nificant higher order corrections 5ji ~ 10%. 



The naive conclusion drawn from such experience |26[ is that one cannot reliably 
calculate the width without ~ 20% uncertainty: 

5r sl /r sl = 6j + 5 n ~ 20% <-> 5\V cb \/\V cb \ ~ 10% . 

On the contrary, theory predicts a strong anticorrelation between 5j and 5jj in a 
consistent perturbative calculation, and that was explicitly checked in [O, 25]. The 



net impact of the calculated (presumably dominant) second-order 0(a 2 ) corrections 
on the value of \Vd,\ appeared to be less than 1% ! Moreover, just neglecting all 
perturbative corrections altogether, both in the semileptonic width and in extracting 
m b from experiment, yields a |V^,| smaller by less than 5% [ 27 . 

Recently, all-order corrections associated with the running of a s in one-loop dia- 
grams (referred to as BLM approximation) were calculated in [p5 ]. Using the most 



accurate model- independent determination of m;, [53], one gets |TB 



W, 



cb\ 



0.0413 



/BR(fi 



0.105 




1.6 ps 



tb 



X 



0.012 



(g - 0.4 GeV 
0.1 GeV 2 




0.006 



5m* b 
30 MeV 



(9) 



The main source of theoretical uncertainty is the exact value of fi 2 (marked with 
a tilde in eq. (H), indicating that a particular field-theoretic definition is assumed), 
which enters through the value of m b — m c , eq. (Q) . A dedicated analysis of the 
lepton spectra will reduce this uncertainty. At the moment a reasonable estimate of 
the uncertainty in ft 2 is about 0.2 GeV 2 , leading to a 2.5% uncertainty in \V c b\. 

The dependence on m b is minor; since we rely here on the well-defined short- 
distance mass ml, there is no intrinsic uncertainty in it. The analysis p3f estimated 
5ml ~ 30 MeV; even for 5ml — 60 MeV, the related uncertainty in | V c b\ is only 1.2%. 

As previously explained, the actual impact of the known perturbative correc- 
tions when relating the semileptonic width to other low-energy observables is very 
moderate, and there is no reason to expect the higher-order effects to be significant. 
With the a priori dominant all-order BLM corrections calculated 25j, one may be 
concerned only with the true two-loop effects 0(a 2 s ). These have not been calculated 
completely yet; however, the recent 0{a 2 s ) calculation [28[ in the small velocity kine- 



matics suggested that they must be small. There are some enhanced higher-order 
non-BLM corrections that are specific to the inclusive widths p9[ . They have been 
accounted for in the analyses M, |16[ , but went beyond those in pO, E5f . Thus it seems 



unlikely that as yet uncalculated second-order corrections can change the width by 
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more than 2-3%; therefore, assigning an additional uncertainty of 2% in \V c b\ is a 
quite conservative estimate. 

Adding up these uncertainties we arrive at 

(S\V cb \/\V cb \)\ th £ 5% . (10) 

The theoretical accuracy in extracting of \V c b\ appears to be better than its cur- 
rent experimental counterpart. This method can be improved further in a model- 
independent way. I think that the 2% level of a defensible theoretical precision can 
be ultimately reached here; an essential improvement beyond that is questionable, 
because of effects of higher-order power corrections and possible violations of duality. 
Similarly, \V u b\ is directly related to the total b —>■ u semileptonic width [fR 



\V ub \ = 0.00458 • [BR(B -> XJv)/0.002] 5 (1.6ps/r B ) 5 . (11) 

Recently, ALEPH announced [[31] a model-independent measurement of the inclusive 
b ^ulv width: BR(B -> XJu) = 0.0016 ± 0.0004. I cannot judge the reliability 
of the quoted error bars in this sophisticated analysis; it certainly will be clarified 
soon. Accepting this input literally, I arrive at the model-independent result 

\V ub \/\V cb \ = 0.098 ±0.013 . (12) 

The theoretical uncertainty in converting T(B — > X u lv) into \V ub \ is a few times 
smaller. 

Let me briefly comment on the literature. It is sometimes stated |2TS|, that the 



uncertainty in T s \ is at least 20%. The origin of such claims is ignoring the subtleties 
related to using the pole mass in the calculations and considering separately the 
perturbative corrections to the pole masses, and to the widths expressed in terms of 
mg*. This is inconsistent on theoretical grounds [pOR , whose relevance was confirmed 
by the concrete numerical evaluations flE| , |25| . The dependence on m b and m b — m n 



used to determine the uncertainty in \ V cb \ was calculated erroneously in (cf. |16|), 
apparently because of an arithmetic mistake that led to a significant overestimate. 
Finally, no argument was given to justify a sevenfold boosting of the theoretical 
uncertainty in m b obtained in the dedicated analysis E3 . 



2.2 Exclusive zero recoil B — > D* tv rate 

Good control of all QCD effects in T s \ was due to the fact that removing constraints 
on the final state to which decay partons can hadronize, makes such a probability 
a short-distance quantity amenable to a direct OPE expansion. A similar approach 
to the exclusive zero-recoil decay rate B — > D* Iv yielded quite an accurate deter- 
mination of \V cb \ as well || |18fl , though with a more significant irreducible model 
dependence and a larger intrinsic uncertainty. The limitation is twofold: constraining 
the decays to a specific final state makes the transition not a genuinely short- distance 
effect; it also suffers from a larger expansion parameter, namely l/m c vs. l/m b . 
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Near zero recoil the decay is governed by the single hadronic formfactor Fp* . In 
the infinite mass limit Fd* = 1 holds; for finite m^ c it acquires corrections: 

Fd* = 1 — (1 — TjA) + 5\/ m 2 + ... . (13) 



The effect of the nonperturbative domain starts with the terms ~ l/(m c , m&) 2 [[33 
|4|, but otherwise is rather arbitrary, depending on the details of the long-distance 
dynamics in the form of wavefunction overlap. This opened the field for speculations 



and controversy f35fl 



The situation as it existed by 1994 was summarized in reviews by Neubert |3| 

r) A = 0.986 ± 0.006 5 1/m 2 = (-2 ± 1)% , (14) 

yielding Fd* — 0.97, and was assigned the status of "one of the most important and, 
certainly, most precise predictions of HQET". Nowadays we believe that the actual 
corrections to the symmetry limit are larger, and the central theoretical value lies 



rather closer to 0.9 || |18[ . Perturbative-wise, it has been pointed out |37[ that the 
improvement |38] of the original one- loop calculation was incorrect, and the proper 
estimate is tja ~ 0.965 ± 0.025; subsequent calculations of the higher-order BLM 
corrections |^5|. |39| confirmed it: tja ~ 0.965 ±0.02. The purely perturbative chapter 
was closed recently with the complete two-loop 0(a 2 s ) result |f28[ r]^ loop ^ = 0.960 ± 



0.007; however, the inherent irreducible uncertainty of the complete perturbative 
series for t/a exceeds the quoted one by a factor of three |25|, |4T], |3| . 

If the mass of the charm quark were a few times larger, in practice the two- 
loop calculation would have been the whole story for Fd*. In reality, the power 
corrections originating from the domain of momenta below ~ 0.6 GeV appear to be 
more significant. Not much can be said about them without model assumptions; 



they have been shown to be negative and exceed about 0.04 p, Q in magnitude. 

The idea of this dynamical approach was to consider the sum over all hadronic 
states in the zero recoil kinematics; such a rate sets an upper bound for the pro- 
duction of D*. This inclusive quantity is of a short-distance nature and can be 
calculated in QCD using the OPE. The result through order 1/m 2 is 



\F D *\ 2 + V l-F-l 2 = Ca(a*) - ^° - ^ n ^ G ( — + — + 2 ) (15) 
* 3m 2 4 ym 2 m 2 3m c m(, J ' 

where Fi are the transition formfactors to charm states i with the mass Mj = Md* + 
ei, and £a is a perturbative factor (the role of /i will be addressed later). Considering 
a similar sum rule for another type of 'weak current', say 02756, yields 

IW=(^-s;)'H--4) («» 



with the tilde referring to the quantities occurring in the transitions induced by this 
hypothetical current. These sum rules (and similar ones at arbitrary momentum 
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transfer), established in || |IBfl , have been subjected to a critical scrutiny for two 
years, but are now accepted and constitute the basis for currently used estimates of 
F D *. 

Since eq. (0) is the sum of certain transition probabilities, it results in a rigorous 
lower bound 

/4 > /i 2 ^ 0.4 GeV 2 . (17) 
The sum rule ([15]) then leads to the model-independent lower bound for S 1 / rn 2: 

r2 Air2 



- 5 1/m 2 > [Mb* - M%) /8m 2 c ~ 0.035 . (18) 

The actual estimate depends essentially on the value of /i 2 . It was suggested in [|J 
to estimate the contribution of the excited states in the l.h.s. of the sum rule (|T5]) 
from to 100% of the power corrections in the right-hand side: 

If so, one arrives at 




-5 1/m 2 = (5.5 ± 1.8)% at /4 = 0.4 GeV 2 

_ S 1/m 2 = (6.8 ± 2.3)% at [l\ = 0.5 GeV 2 (20) 

-5 1/m 2 = (8.1 ± 2.7)% at /i 2 = 0.6 GeV 2 

The QM meaning of the sum rules is transparent The act of a semileptonic 
decay of the b quark is its instantaneous replacement by a c quark. In ordinary QM 
the overall probability of the produced state to hadronize to something is exactly 
unity. Why are there nonperturbative corrections in the sum rule? The answer is 
that the 'normalization' of the weak current 07^756 is not exactly unity and depends, 
in particular, on the external gluon field. Expressing the QCD current in terms of 
the nonrelativistic fields used in QM one has, for example, through order 1/m 2 : 

C7fc75& <-> crfc - ( Tr-^{0iD) 2 a k + -^a k (aiD) 2 - — ^ — (aiD)a k (aiD) ) . (21) 
\Qm^ 8m£ 4m c mi) J 

The last term just yields the correction seen in the r.h.s. of the sum rule. Let me note 
that in the standard HQET analysis, the first two terms in the brackets are missing 
(see, e.g., Ref. |35| ) and the dominant effect ~ l/"^ 2 is lost; the nonrelativistic 
expansion was correctly done in the works by the Mainz group [43|. 



The inequality fi 2 > \i G in QM expresses the positivity of the Pauli Hamiltonian 
^(aiD) 2 = ^((iD) 2 — |crG) [|]. It is interpreted as the Landau precession of a 
charged (colored) particle in the (chromo)magnetic field where one has (p 2 ) > \B\. 
Although the QM average of B in the B meson is suppressed, the chromomagnetic 
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field is proportional to the spin of the light degrees of freedom and is thus essentially 
non-classical, which enhances the bound and makes up for the suppression. 

The perturbative factor £4 is not equal to rj\ |T8], ^(J and depends on the sep- 
aration scale /1. Unlike T)a, which in principle cannot be defined theoretically with 
better than a few percent accuracy, £a(aO i s well-defined; no significant uncertainty 
is associated with it, £4 ~ (0.99) 2 . 

Allowing a very moderate variation of // 2 between 0.4 GeV 2 and 0.6 GeV 2 only, 
we see that — 5 x / m 2 varies between 3.5% and 11%; moreover, since there are no 
model-independent arguments to prefer any part of the interval, the whole range 
must be considered equally possible. Adding small perturbative corrections we end 
up with the reasonable estimate Fd* ~ 0.9 . It is curious to note that at a 'central' 
value x — 0-5 the dependence of the zero-recoil decay rate on /i 2 through 8y m 2 effect 
practically coincides with that of T a \{B) (see eq. (||) ) although they actually vary 
in opposite directions. The typical size of the l/mn corrections to the exclusive 
zero-recoil decay rate is thus significant, around 15%, which is expected since they 
are driven by the scale m c ~ 1.3 GeV. It is evident that \jw? c corrections in Fd* not 
addressed so far are at least about |(0.15) 3 / 2 ~ 2-3%.Q 

Thus, I believe that the current theoretical technologies do not allow to reliably 
predict the zero recoil formfactor Fd* with a precision better than 5-7% in a model- 
independent way; its value is expected to be approximately 0.9, although a correction 
to the symmetry limit twice smaller, as well as larger deviations, are possible. It 
is encouraging that the 'educated guess' Fd* ~ 0.9, which emerged from the first - 
and so far the only - dynamical QCD-based consideration || [LB], yielded a value of 
\V c b\ close to a less uncertain result obtained from T S \(B). 

Future, more accurate data will enable us to measure Fd* with a theoretically 
informative precision using \V c b\ from T S \(B), and thus provide us with deeper insights 
into the dynamics of strong forces in the heavy quark system. 

Certain statements in the literature deserve comments. Neubert claimed []4"0"|, 
41] that the sum rule (|l|) cannot be correct, since 1/m 2 renormalons allegedly 
mismatch in it. It was failed to realized in these papers that in Wilson's OPE the 
IR renormalons are always absent from any particular term. On the other hand, IR 
renormalon calculus can still be applied if the OPE relation is considered in the pure 
perturbation theory itself, and formally setting /1 = 0. In particular, this amounts 
to subtracting a 'perturbative piece' from the observable probabilities. However, in 
this way the perturbative terms appear in the left-hand side as well, and these terms 



were ignored in Refs. [40, 41 



It was suggested in Ref. [[H] that the higher-order radiative corrections to the 
sum rules are too large and allegedly make them next to useless. Such conclusions 
emerged from a theoretically inappropriate treatment. The concrete numerical anal- 
ysis in the OPE quoted below, on the contrary, suggests a quite moderate impact 



1 This is consistent with the fact that the 1/rriQ IR renormalon ambiguity in r) 2 A constitutes 5% 



at Aqq D ~ 220 MeV @ 
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of radiative corrections. According to []4"4" |, the perturbative corrections to the sum 
rule of the type of eq. (|T6| ) weaken the bound for the expectation value of the kinetic 
operator to such an extent that it becomes non-informative. One must realize that, 
in the quantum field theory, the renormalized operators can be defined in different 
non-equivalent ways; — Ai addressed in [ii] is known to be different from fi 2 . More- 
over, the only field-theoretic definition of the kinetic operator Q(iD ) 2 Q given so far 
was made in [JTSJ and, for it, fi 2 > [i G always holds. As for — Ai, a parameter in 
HQET, its definition beyond the classical level has never been given; the procedure 
adopted in Ref. ||4j reduces to an attempt to completely subtract the 'perturbative 
piece' of 



2 ,,2 



Ai = //^(/i) - ci (a s (//)/7r)/i - c 2 (a s (//)/7r) fi 



(22) 



(the method to calculate q was elaborated in [JTJJ). Yet it has been known for a 
long time that such a program theoretically cannot be performed: the series in 
eq. (|2"2"D is factorially divergent and cannot be assigned a meaningful number. No 



wonder the second-order BLM correction calculated in [44| seemed to be dangerously 
large: there can be no bound established for a quantity that is not defined. Moreover, 
the situation is clear in the BLM approximation, where all q are readily calculated: 



the series, whose second term was discussed in |]44 |, is divergent and sign- varying, 
so using merely the second term is misleading for any numerical estimate. 

The above subtleties are peculiar to the field-theory analysis. Inequality fil > \j? G 
must hold in any QM model relying on a potential description without additional 
degrees of freedom, if the heavy quark Hamiltonian is consistent with QCD. Unfor- 
tunately, a failure to realize this fact is seen in a number of recent analyses. 

A second-order BLM analysis of the sum rule ( |15"D for Fd* was also attempted 



m 



||4]] and claimed to destroy its predictive power (the first-order calculation had 
been performed in [113]). This calculation, however, was not done consistently, and 



the actual effect is smaller f42fl . Let me define J7a(a*) = ^a(^) 1 ^ 2 ] the quantity ?7a(/x) 
must be added to <5i/ m 2 in the framework of the OPE instead of r] A in the model 



calculations. Then, at a reasonable choice [i ~ 0.5 GeV, Aq2 d = 300 MeV, one has 



Va(h) = 1 tree level 

Va(i^) = 0.975 one loop 
rj A (p) = 0.99 all-order BLM 



(23) 
(24) 
(25) 



Clearly, the effect of the calculated perturbative corrections is not drastic and 77a(aO 
is very close to the value of 0.98 adopted in the original analysis [§. 

Conceptually, the deficiency of the alternative application of the original sum 
rules of Refs. || |18[ adopted in [(44| , is to gauge £4 on the value of rj\ as it has been 
defined in the HQET (the idea that £4 is to be identified with rj\ ascends to [40"|). 
However, it is tja that is ill-defined, and only for this reason must the difference 
between the stable Wilson coefficient £4 and r)\ suffer from large corrections. It is 
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worth noting that, in reality, r/ a cannot be equal to a matching coefficient of 07^756 
to a corresponding current in any effective field theory. 

Smaller theoretical uncertainty, 2.5% and 3%, is now quoted by Neubert for Su m 2 



and Fd*-, respectively. The former was obtained in Ref. ||40|| , in what he calls a "hy- 
brid approach", which reduces to assigning the fixed value ft 2 = 0.4 GeV 2 and using 
it in the sum rule ( |I~5[) within the same model assumption of eq. ( |l~9"D as suggested in 
Hl : < x < 1-Q Correspondingly, the quoted number for d>i/ m a practically coincided 



with the first line of eqs. (|20|). In reality, allowing to vary within any reason- 



able interval significantly stretches the uncertainty. Moreover, the analysis [40 



was 



based on using t]a as a perturbative factor assuming, literally, that in the proper 
treatment the final result would not be changed numerically - which is just the case 
according to Ref. On top of that, the uncertainty in the definition of t\a due to 
1/m 2 and 1/m 3 IR renormalons constitutes 2-3% each and is an additional one in 
the adopted usage of the sum rules. Altogether, the stated theoretical confidence of 
those estimates cannot be accepted as realistic. 

Recently, Ref. |^5j claimed to have established an intriguing relation between the 
slope and the curvature of the formfactor near the zero recoil point, using analyticity 
and unitarity of the amplitudes. If correct, this would reduce the experimental 
uncertainties in extrapolating the rate to the zero recoil. However, both the below- 
threshold contribution to the dispersion integral and the l/m c power corrections 
4Tf| to the heavy quark symmetry relations were grossly underestimated; therefore 



the relation stated in [45| rather should not be used for deriving model-independent 
experimental results. 



3 The semileptonic branching fraction 

The QCD-based HQE provides a systematic framework for calculating the total 
widths of heavy flavors, which are not amenable to the traditional methods of HQET. 
The difference between nonleptonic and semileptonic widths appears only at a quan- 
titative level. The only assumption is that the mass of a decaying quark (actually, 
the energy release) is sufficiently large; for a review, see Wft . 

The overall semileptonic branching ratio BR s i(S) seems to be of a particular 
practical interest: while the simple-minded parton estimates yield BR s i(S) ~ 15% 
]48| , experiments give smaller values BR s i(5) ~ 10.5-11.5%. The leading 1/tuq 



effects in the nonleptonic widths were calculated in |7|, 0, 13| ; they exhibit some 



cancellations and one literally gets |49j a downward shift ^ 0.5%. Estimated 1/m; 



corrections do not produce a significant effect either f)0| |. As a result, most of the 



attention was paid to a more accurate treatment of the perturbative corrections, 



2 I disagree with the statements of Q, reiterated in later papers, suggesting that the original 
analysis ]8|, |l8) missed some elements of the heavy quark spin- flavor symmetry; on the contrary, 
it was stated in the latter paper that all these relations automatically emerge from the sum rules 
that replace the QM wavefunction description in the quantum field theory. 
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including the effect of the charm mass in the final state. It was found |51| , |52|| that 
the nonleptonic width is indeed boosted up. 

Since the inclusive widths are expanded in inverse powers of energy release, one 
expects larger corrections or even a breakdown of the expansion and violation of 
duality in the channel b — > ccs(d); however, this channel can be isolated via charm 
counting [[/J. The original experimental estimate n c ^ 1.15 did not allow one to 
attribute the apparent discrepancy to it, and gave rise to the so-called 'BR sl versus 
ric problem. 

The perturbative corrections in the b — > cud itself cannot naturally drive BR s i 
below 12.5%; the calculation for the b — > ccs(d) channel is less certain and, in 
principle, admits increasing the width by a factor of 1.5-2, leading to n c ~ 1.25-1.3. 
In the latter value of BR s i as low as 11.5% can be accommodated. 

The experimental situation with n c does not seem to be quite settled yet: n c = 
1.134 ± 0.043 (CLEO), n c = 1.23 ± 0.07 (ALEPH). In a recent analysis |53| it was 
argued that consistency requires a major portion of the final states in b — > ccs to 
appear as modes with kaons but D s , which previously escaped proper attention. 
This allows for a larger value n c ~ 1.3 needed to resolve the problem with BR s i. 
The dedicated theoretical analysis |54| shows that, indeed, the dominance of such 
modes is natural and does not require violation of duality. Thus, if a larger value of 
n c ~ 1.25 is confirmed experimentally, the problem of BR s i will not remain. 

In my opinion, however, we cannot consider even this successful scenario as a 
complete QCD-based theoretical prediction of BR s i; a strong enhancement of a tree- 
level-unsuppressed channel raises doubts about the trustworthiness of its one-loop 
calculations. The possibility to get the necessary enhancement should be rather 
viewed as an indication of the presence of large effects working in the right direction. 



4 Lifetimes of beauty particles 

A thoughtful application of the HQE to charm lifetimes demonstrated that the 
actual expansion parameter appeared to be too low to ensure a trustworthy accurate 
description, so that a priori one expects only emergence of the qualitative features. 
Surprisingly, in most cases the expansion works well enough even numerically (for a 



recent review, see ||55]| ). 

Applying the expansion to beauty particles one expects a decent numerical ac- 
curacy, although the overall scale of the effects is predicted to be small, making a 
challenge to experiment: 



t b -/t b o 
Tb s /t b o 
(t b l ~ t b s)/t Bs 
t aJt b o 



Sb 



[50] 



1 + 0.04 ^ 180MeV 
1+0(1%) 150 
0.18 
0.9 



200 MeV 



m 



EXP: 
EXP: 

EXP: 



1.04 + 0.04 
0.97 + 0.05 

0.78 + 0.06 



(26) 
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These differences appear mainly as \jm\ corrections and, depending on certain four- 
fermion matrix elements, cannot be predicted at present very accurately, in partic- 
ular in baryons. For mesons the estimates are based on the vacuum saturation 
approximation, which cannot be exact either. The impact of non-factorizable terms 
has been studied a few years ago in |57J and possibilities to directly measure the 
matrix elements in future experiments were suggested. 

The apparent agreement with experiment is obscured by reported lower values of 
T\ b . Since the baryonic matrix elements are rather uncertain, a few model estimates 
have been done [p8 |. All seem to fall short; however, this might be attributed to 
deficiencies of the simple quark model. Nevertheless, it was shown |59[ that irre- 
spective of the details one cannot have an effect exceeding 10-12% while residing in 
the domain of validity of the standard l/rriq expansion itself; the natural 'maximal' 
effects that can be accommodated are ~ 7% and ~ 3% for weak scattering and 
interference, respectively. 

Thus, if the low experimental value of T\ b is confirmed, it will require a certain 
revision of the standard picture of the heavy hadrons and of convergence of the 1/rriQ 
expansion for nonleptonic widths with the offset of duality in beauty particles. 

Recently, the problem of the accuracy of the calculations of 5t b based on the 
factorization was emphasized again in [EjJ. It is difficult to agree, however, with 
the wide intervals, up to ±20% allowed for the difference between t b + and r B o; the 

1571 



constraints discussed in 



were missed. One can see that the values of hadronic 



parameters saturating such large differences would move one beyond the domain of 
applicability of the whole expansion used in [BO . 



5 1/?tiq expansion and duality violation 

Duality violation attracts more and more attention in the context of the heavy 
quark theory; a recent extensive discussion was given in [pjj. The expansion in 
1/tbq is asymptotic. There are basically two questions one can ask here: what 
is the onset of duality, i.e. when does the expansion start to work? The most 



straightforward approach was first undertaken in |62|| , and no apparent indication 
toward an increased energy scale was found. Another question, of how is the equality 
of the QCD parton-based predictions with the actual decay rates achieved, was rarely 
addressed. An example of such a problem is easy to give. 

The OPE states that no terms ~ 1/^q can be in the widths and the leading 
terms start with l/roq. However, the OPE per se cannot forbid a scenario where, 
for instance, 

ST Hq /T Hq ~ C sin (m Q p)/(m Q p) , p ~ Aq^ d . (27) 

In the actual strong interaction, m,b and m c are fixed, so from the practical viewpoint 
these types of corrections are not too different - but the difference is profound theory- 
wise! It is a specific feature of the OPE in Minkowski space, and it can hardly be 
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addressed, for example, in lattice calculations. Their complete control requires a 
deeper understanding of the underlying QCD dynamics beyond the knowledge of 
first few nonperturbative condensates. 



The literal corrections of the type of eq. (27) are hardly possible; the power 
of 1/toq in realistic scenarios is larger, and they must be eventually exponentially 
suppressed though, probably, starting at a higher scale f6~lfl . But a theory of such 



effects is still in its embryonic stage and needs an experimental input as well. 

The possibility has been discussed for some time |x| to have certain unidentified 
corrections to the (nonleptonic) widths, eventually leading to the dependence 

T n 4 ~ M\ Q (28) 

and manifesting an explicit 1/rriQ effect: 

6T Hq /T Hq ~ 7 SM Hq /M Hq , 7^5. (29) 

Such scaling in the intermediate energy domain cannot literally contradict OPE if 
the offset of duality has not been passed yet. But is this possibility natural? 

Leaving aside the QCD-based arguments completely, one must still account for 
the charm mass in the final state, and thus differentiate between Mp and Ma c in 
the decays of B and A&, respectively. Although M^ b is notably larger than M#, Ma c 
exceeds by almost the same amount! Counting only the phase space factors in 
analogy with the free quark decay one would get for the b — > c transitions 7^2 
which thus seems to be a more natural value in the 'poor man on the street' 



hypothesis |]] considered in |]63|] . It is worth noting that the fit of 7 in charm par- 
ticles does not convincingly indicate favoring 7 = 5: incorporating the calculated 
and \jrn? c corrections destroys it - while discarding them as a part of the 'poor 
man' philosophy makes it impossible to explain the very different values of td and 
BR s i(D). The short A& lifetime thus still seems to constitute an important problem 
for the whole heavy quark theoretical community. 
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